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Abstract 



Predictions for Deep Virtual Compton Scattering are obtained in a two com- 
ponent dipole model of diffraction. The model automatically includes hard and 
soft components and implicitly allows for "hadronic" contributions via large 
dipoles. It is also applicable to real Compton Scattering, which provides an 
important constraint. 
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1 Introduction 



Deep Virtual Compton Scattering (DVCS) on protons, j*p — > jp, is an im- 
portant reaction for the study of diffraction. In the standard perturbative 
QCD approach the amplitude is described by skewed parton distributions jl| 
for Q 2 > Ql, for some non-zero Qq. The skewed parton distributions corre- 
spond to operator products evaluated between protons of unequal momenta. 
These are generalisations of the familiar parton distributions of deep inelas- 
tic scattering, and like them satisfy perturbative evolution equations |g{ which 
enable them to be evaluated at all Q 2 in terms of an assumed input at some 
appropriate Q 2 — Q 2 . 

Preliminary data have been presented || which are consistent with QCD predic- 
tions H . The QCD calculations require the input skewed parton distributions 
at the reference point Qq. In 0] these are obtained by estimating their ratio 
to "ordinary" parton distributions at Qq = 2.5 GeV 2 using arguments based 
on the aligned jet model Js|. In practice this is equivalent |6) to the simplest 
diagonal generalised vector meson dominance model 0. The role of vector 
mesons, particularly the p, in providing a "hadronic" contribution to DVCS via 
the vector-meson-dominance mechanism ^*p — > Vp, V — > 7 has been calculated 
in |6| and shown to be important at values of Q 2 well beyond the input reference 
point Qq. 

The model of p electroproduction used in Q is a two-component model, combin- 
ing soft and hard contributions which are respectively primarily non-perturbative 
and perturbative in origin. There is considerable evidence in 7*73 interactions 
that the nominally perturbative regime can be strongly influenced by non- 
perturbative effects. This is an obvious feature of dipole models of deep inelastic 
scattering || — [[li) , where the contribution from large dipoles can extend to sig- 
nificantly large values of Q 2 ; and in two-component models [^5) — 0|, in which 
there is a non-perturbative term at all Q 2 by construction. 

Here we present the results of a model which incorporates these effects not in 
the language of parton distributions but rather that of dipole cross sections. In 
the framework of functional integration a two-component dipole cross section is 
constructed based on the approach of Jll], [l8|. The model covers the complete 
range from the real 7^ cross section through DVCS, with one photon virtual, to 
deep inelastic scattering, with both photons virtual. The model as applied to 
DVCS contains no free parameters at t — t min : they have already been deter- 
mined from the pp and 7*p total cross sections. To obtain the integrated cross 
section requires knowledge of the logarithmic slope of the differential cross sec- 
tion. This can also be determined from the model, and for p electroproduction 
[ pp| is in good agreement with the experimental value over the relevant range 
of Q 2 . We use this as a reasonable estimate in the present dipole approach. 
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2 The model 



We present a calculation for the cross section of Compton scattering with differ- 
ent virtualities of the photons. The treatment is based on a functional approach 
of high energy scattering with small momentum transfer |2lj| in which the scat- 
tering of a qq pair is expressed as the functional average of a Wegner- Wilson 
loop with light-like sides ^3j . The functional average is performed in a 
model for non-perturbative QCD [25). The model depends essentially on 
two typically non-perturbative parameters: the strength of the gluon correla- 
tor, given through the gluon condensate (g 2 FF), and its correlation length a. 
The linear confining potential is also defined in terms of these two parameters 
p5| . Nucleons are most simply treated as a quark-diquark system, that is 
effectively as a dipole, although this is not essential. The two parameters of 
the model were obtained by fitting the isoscalar part of pp and pp scattering at 
W = 20 GeV, and were found to be (g 2 FF) = 2.49 GeV 4 and a = 0.346 fm. 
These values are within the range determined from lattice calculations, and give 
the correct value for the slope of the confining potential. 

In this approach scattering amplitudes with small momentum transfer can also 
be calculated, but in this note we concentrate on forward scattering. Then t he 
approach becomes a simple dipole model and the Compton scattering amplitude 
can be expressed as an integral over the product of the dipole-proton cross 
section ad(R) and the overlap of the photon wave functions with virtualities Q\ 
and Q\ and helicity A, pt{Q\, Q\, R, z): 



T(s, t = 0) = 2tt J dz J ' dRRp*(Ql,Ql,R, 



z)a d {R) (1) 



where z is the longitudinal momentum fraction of the quark in the photons and 
R is the radius of the quark-antiquark dipole, the weak dependence of the dipole 
cross section Q is neglected. The dipole cross section ad{R) as evaluated in 
the model can be very well described by 

a d (R) = 0.098(( 5 2 J FF)a 4 ) 2 i?(l - exp ( - — )) (2) 

The dimensionless constant (g 2 FF)a 4: has the numerical value 23.8. If a and R 
are measured in fermi, the result is in millibarn. 

As it stands the model has no energy dependence. The increase with energy of 
hadronic cross sections as W 2esott , where the intercept of the pomeron trajectory 
is I + e so f t with e so f t ~ 0.09 ± 0.01, can be incorporated in two ways. Either 
one lets the radii of hadrons increase with W 2 (23) , [^6| — , or one takes the 
model as a determination of the coefficient of pomeron exchange and includes a 
factor (W/W ) 2€ °° {t with W = 20 GeV. These two approaches give very similar 
results and we adopt the latter in this paper as it is the more convenient in the 
present context. 
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To incorporate the greater energy dependence observed in reactions with a hard 
scale, such as deep inelastic scattering, the two-pomeron approach of |1| was 
adapted to the model in jl6j and applied successfully to the photo- and electro- 
production of vector mesons and to the proton structure function over a wide 
range of Q 2 . As in |L5| it was found that the soft-pomeron contribution to 
the proton structure function initially increases with increasing Q 2 , has a broad 
maximum in the region of 5 GeV 2 and then decreases slowly as Q 2 decreases fur- 
ther. In the context of the present model this is a consequence of the decreasing 
interaction strength with decreasing dipole size. 

It is worth recalling the salient features of this version of the two-pomeron model 
to illustrate the distinction between the soft- and hard-pomeron contributions. 
Following |l6| it is assumed that all dipole amplitudes in which both dipoles 
are larger than the correlation length a = 0.346 fm are dominated by the soft 
pomeron, and the energy dependence is given by (W/Wo) 2e ° olt . This ensures 
that the hard pomeron has essentially no impact on purely hadronic scattering. 
If at least one of the dipoles is smaller than a fm then the energy dependence 
is replaced by (VF/Wo) 2ehard , with e h ard = 0.42. This is the value found by 
[ [l5| and is larger than the value used in pl| . This difference in ehard = 0.42 
between the present paper and |^6) is accounted for by a different treatment 
of the z integration. In Jl6| there was a cutoff in the z integration, requiring 
that W z (1 — z) > 0.2 GeV in order to ensure that the quarks in the centre- 
of-mass frame have at least this energy. This z cutoff introduces an additional 
energy dependence, particularly for highly virtual photons. In our approach we 
follow the general lines of the dipole models and integrate z from zero to one 
and therefore need a higher intercept for the hard component. It was found 
in [jl6| that the extrapolation of the non-perturbative model to high values of 
Q 2 overestimates the small dipole contributions and therefore the dipole cross 
section is put to zero below a radius R c = 0.16 fm. 

The quark-antiquark overlap density of photons with virtuality Q\ and Q\ re- 
spectively and helicity A in lowest order perturbation theory is 

IN a 

p°{Qj,Q 2 ,R,z) = —^-e 2 f Q 1 Q 2 z 2 (l-z) 2 K (e 1 R)K (e 2 R) 

p^(,Ql,Ql,R, z) = Q ^e 2 f ({z 2 + {l-z) 2 )€ 1 K 1 (e 1 R)e 2 K 1 {e 2 R) 

+m 2 K (eiR)K (e 2 R)y (3) 

Here e/ is the charge of the quark in units of the elementary charge; my is the 

Lagrangian quark mass; = \J z{\ — z)Q 2 + m 2 ; A = indicates a longitudinal 

photon and A = ±1 indicates a transverse photon. The singularities of the 
Bessel functions at R = do not cause any problems, since for the evaluation 
of observable amplitudes the density is multiplied by the dipole cross section 
which is proportional to R 2 for small values of R. 

For longitudinal photons the factor z 2 (\ — z) 2 in the density p* ensures that 
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Figure 1: The integrated cross section for the reaction j* p — > jp as function of 
the virtuality Q 2 of the incoming photon at an average (W) — 75 GeV 



the main contribution comes from the region z « |, and the relevant scale 
is \Q 2 + m 2 . Thus contributions from large dipoles are suppressed at large 
Q 2 . However for transverse photons the endpoints z — and z = 1 are not 
suppressed and large dipoles can contribute to the perturbative regime even if 
the virtuality is quite high. The z cutoff of |Q mentioned above is therefore 
most important for large values of Q 2 . 

The results (g) are reliable for large values of Q 2 and/or large quark masses. 
For small values of Q 2 and light quarks the separation of the quark-antiquark 
pair can become large and confinement effects become important. An effective 
way in which to deal with this problem is to introduce a Q 2 -dependent effective 
constituent quark mass as an infrared regulator in the photon wave function. 
This procedure is justified in p9{ where the value of the light quark effective 
mass at Q 2 — is determined from the two-point function of the vector current 
as iriQ q = 0.21 ± 0.015 GeV and for the strange quark as mg s « 0.31 GeV. Since 
in this note we consider only processes where one photon is real we always use 
these effective mass values at Q 2 = 0. The charmed quark mass was taken as 
m c = 1.3 GeV. 



3 Results 



A convenient fit to the model calculation for < Q 2 < 100 GeV 2 (good in the 
amplitudes up to 1 % for Q 2 < 20 GeV 2 and 3 % for higher Q 2 values) for one 
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Figure 2: The integrated cross section for the reaction j* p — > jp as function of 
the center-of-mass energy W at an average virtuality of the incoming photon of 
(Q 2 ) = 4.5 GeV 2 



photon real, the other one with virtuality Q 2 is given by: 

— (Q 2 ,W,t = t min ) = 
at 



l^{r soh (Q 2 )(W/W ) aW + n iatd (Q 2 )(W/W ) 



0.84 
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where W = 20 GeV and 



311 /zb/GeV 2 (4) 



rWt(Q') = (5.33 - 1.33 exp(-4Q 2 /Q 2 ) + 5.37 Q 2 /Q 2 
r h ard(Q 2 ) = (49.42 - 7.65 exp(-4Q 2 /Q 2 ) + 4.94 Q 2 /Q 



(5) 



with Qq = 1 GeV 2 The present calculation provides the j*p — > jp amplitude 



only for forward scattering, that is at t — t„ 



proton 



- /s , so it is nec- 



essary to make an assumption about the i-dependence to obtain an integrated 
cross section. As the process is dominated by light-quark dipoles, it is reason- 
able to take the i-dependence to be similar to that for p electroproduction. For 
defmiteness we take the logarithmic slope to be B — 7 GeV -2 , which is the 
average value over the Q 2 range of the preliminary HI data ||. Such a value of 
the average slope is also the result of model calculation for electroproduction of 



p mesons Ed, 31 



The relation between the ep — > epj cross section and the j*p — > "fp cross section 
is M 



dWdQ 2 



W 



7T Q 2 (W 2 + Q 2 



proton 



(1 + (1 - yf)a Tp ^ 



7P) 



(6) 
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with y = (W 2 + Q 2 — m 2 Ioton )/(s — m 2 Ioton ). Here y/s is the centre-of-mass 
energy of the ep system and W is that of the j*p system. 

The Q 2 dependence and energy dependence of the predicted integrated j*p —> 
"fp cross sections are given in figures [l] and [2] respectively. The preliminary HI 
data H are for ep — > ep7 and contain a large background from the QED Bcthc- 
Heitler process. However as the latter is purely real and as the QCD amplitude 
is mainly imaginary (the real part is expected to be only about 10 — 15% of the 
imaginary part), interference is not large and the Bethe- Heftier cross section 
can simply be subtracted. A comparison with the preliminary HI data after 
subtraction of the Bethe-Heitler is made in figures |l| and |[ where we have 
converted the experimental ep — > epj to 7*p — > jp, after subtraction of the 
Bethe-Heitler cross section, using (^). The only serious discrepancy between 
the model and the preliminary data is at Q 2 = 3.5 GeV 2 . This is reflected in 
the normalisation of the integrated cross section in figure || as the low-Q 2 point 
dominates this cross section. We simply note that the inclusion of the real jp 
data is seen to provide an important constraint on models. 
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